Anomalous Scaling of Fracture Surfaces 
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We argue that fracture surfaces may exhibit anomalous dynamic scaling properties akin to what 
occurs in some models of kinetic roughening. We determine the complete scaling behavior of the 
local fluctuations of a brittle fracture in a granite block from experimental data. We obtain a global 
roughness exponent \ = 1-2 which differs from the local one, xioc = 0.79. Implications on fracture 
physics are discussed. 
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The study of the morphology of fracture surfaces is 
nowadays a very active field of research. From the early 
work of Mandelbrot et. al. pL much effort has been paid 
to the statistical characterization of the resulting fractal 
surfaces in fracture processes. Scale invariance has been 
found in many experiments Jl]-|j| and it is now well estab- 
lished that, in general, crack surfaces exhibit self-affmc 
scaling properties in a large range of length scales (see 
JlO| for a more detailed account of experiments). 

A self-affmc surface /i(x) is invariant under an 
anisotropic scale transformation, in the sense that /i(x) 
has the same statistical properties as A~ x ft,(Ax), where \ 
is the roughness exponent. Initially, the hope in studying 
the surface morphology was to relate geometry to me- 
chanical properties (toughness, plasticity, etc) in order to 
obtain a material characterization by means of the rough- 
ness exponent. However, experimental results in very dif- 
ferent types of materials (from ductile aluminium alloys 
to brittle materials like rock) seem to support the idea 
of a universal roughness exponent which is very indepen- 
dent of the material properties. For three-dimensional 
fractures, an exponent x(3D) ~ 0.8 — 0.9 [jl]-|| has been 
measured, whereas in dimension two x(2D) ~ 0.6 — 0.7 
[^-^|. It seems reasonable to expect that material prop- 
erties should affect the fracture surface roughness. In 
particular, toughness and anisotropy should be relevant 
for the fracture crack morphology. However, the above 
mentioned experimental results seem to lead to the sur- 
prising conclusion that there is no correlation of \ with 
mechanical properties. 

The treatment of the fracture crack as a self-affinc 
rough surface leads in a natural way to the close field 
of kinetic roughening (see jll[| for recent reviews in the 
subject). A direct mapping of the crack at the stationary 
state into the Kardar-Parisi-Zhang equation has been 
proposed in Ref. |X3f| . The crack surface has also been 
considered as the trace of the crack front whose propa- 
gation is modeled by other types of non-linear Langevin 
equations Q. 

In this Letter, we will show that crack surfaces have 
very much in common with those obtained in growth pro- 
cesses exhibiting what is called an anomalous dynamic 
scaling. We argue that the scaling of the fluctuations of 
crack surfaces is, in the sense of kinetic roughening, in- 



trinsically anomalous rather than a simple Family- Vicsek 
one. The analysis of the roughness data of a crack exper- 
iment demonstrates the validity of our argument. Phys- 
ical consequences of this scaling for fracture will be dis- 
cussed. 

Family- Vicsek scaling. -In practice, the self-affine char- 
acter of a surface in dimension d + 1 is shown by study- 
ing the scaling of the fluctuations of the surface height 
over the whole system of total size L. The invariance 
property under the scale transformation implies that the 
global width at time t, W(L,t) = {(h(x,t) ~h{t)) 2 ) 
where x is the position, has to scale as p~5| 
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W(L,t) = L*f(L/t 1 /*), 



(1) 



where z is the dynamical exponent and / (it) is the scaling 
function 



,, s i const, if m <C 1 
JW ~ i if u > i. 



(2) 



X is the roughness exponent and gives the scaling of the 
surface in saturation, W(L, t 3> L z ) ~ L x . These two 
exponents characterize the universality class of the par- 
ticular growth model. Equivalently, the scaling behavior 
of the surface might be obtained by looking at the local 
width over a window of size I <C L: 



•j(l,t) 



t x ' z if t < r 

l x if t > P. 



(3) 



This is the method actually used in real experiments 
where the size of the system L remains constant and the 
fluctuations are calculated over scales I <C L. Note that 
the local width saturates at time I* as w(l,t ^S> l z ) ~ l x 
independently of the system size. 

A complementary technique to determine the critical 
exponents of a growing surface is to study the Fourier 
transform of the interface height in a system of lin- 
ear size L, h(k,t) — L~ d / 2 X)x[M x > *) — h(t)] exp(ikx), 
where the spatial average of the height has been sub- 
tracted. In this representation, the properties of the sur- 
face can be investigated by calculating the power spec- 
trum S(k, t) = (h(k, k, t)}, which contains the same 
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information on the system as the local width. In most 
growth models, the power spectrum scales as 



S(k,t) = k-Px+Qaikt 1 /'), 
where s is the simple scaling function 



s(u) 



const, if u 3> 1 
u 2x+d if u < 1. 



(4) 



(5) 



This form for the power spectrum can be easily inverted 
to obtain the scaling behavior of both local and global 
widths described above. 

Anomalous scaling and implications.- In very recent 
studies of several growth models , it has been found 

that the surface fluctuations may exhibit an anomalous 
scaling, in the sense that, although the global width be- 
haves as in (0)-(||) local surface fluctuations do not satisfy 
Eq.(pl), but scale as 



w(l,t) 



t x ' z if t < 

f/3, IXioc if t ^ 



(6) 



where the exponent /3* = (x — Xioc)/ z is an anomalous 
time exponent and xioc the local roughness exponent. 
Thus, in the case of anomalous scaling two exponents, 
Xioc and X: enter the scaling and must be taken into 
account to give a complete description of the scaling be- 
havior of the surface. An outstanding consequence is that 
the local width does not saturate at times l z but when 
the whole system does, i.e. at times L z , giving an un- 
conventional dependence of the stationary local width on 
the system size as 



w(l,t > L z ) ~ l^"c L x- 



(7) 



in such a way that the magnitude of the roughness over 
regions of same size I at saturation is not just a func- 
tion of the window size but also of the system size, which 
is distinctly different from what happens in the Family- 
Vicsek case. 

Owing to several experimental limitations, anomalous 
scaling is difficult to observe (see for kinetic roughen- 
ing experiments in which anomalous scaling was found). 
Since the system size L of experiments can hardly be 
changed over a broad range, the dependence of the global 
width on the system size cannot be actually determined. 
Only local fluctuations over a window I can be measured. 
Moreover, very often the time evolution of a crack can- 
not be monitored and in most experiments only the final 
crack surface is analyzed, i.e. Eq.(Q) for a fixed system 
size L. This immediately leads to the conclusion that, 
whether anomalous scaling exists, only xioc is actually at 
reach of the methods currently used in experiments. 

In terms of the power spectrum, the existence of a lo- 
cal exponent xioc ^ X comes from a nonstandard form of 
the scaling function s(u) in Eq.([|). It has recently been 
shown [|l9| that the anomalous scaling (JsJ) is associated 
with either super-roughening EQ] or a power spectrum 



that satisfies the dynamic scaling behavior stated in (Q) 
but with a distinct scaling function 

f u 2 <X-Xloo) if u > i 

S{U) ~\ u 2x + d ifu«l. (8) 

So, in the stationary regime (at times t ^S> L z ) the power 
spectrum scales as S(k,t) ~ k~( 2x,oc+d ^ L 2 ( x ~ Xloa \ and 
not simply as k~( 2x+d ^ as corresponds to a standard scal- 
ing. This means that experimental determinations of the 
roughness exponent of fracture crack surfaces from the 
decay of the power spectrum with fc also give a measure 
of xioc and not x- 

Most of the experimental studies are unable to follow 
the crack in time and much important information about 
the complete scaling is lost. In the majority of the ex- 
periments one has to deal with a static fracture surface 
and fluctuations of its height are evaluated over windows 
of different sizes I. So, neither ^) nor ^) scaling forms 
are actually tested in fracture experiments. 

Experiment. - In the following we present an analysis of 
the data describing growth of the crack roughness. In this 
experiment a fracture was initiated from a straight notch 
in a granite sample (25cmx25cmx 12cm) ||. It is a mode 
I unstable crack. The crack roughness increases from two 
hundredth of a millimeter to several millimeters. To- 
pographies of two areas of 5cm x 4cm were recorded with 
a first mechanical profiler along a regular grid (100 paral- 
lel profiles). The x direction which is parallel to the initial 
notch was sampled with 1050 points (Ax = 50/im). The 
grid step along the perpendicular direction (i.e. the crack 
propagation direction) was Ay = 350/im. A third map 
(5cm x 5cm) was obtained from a second and independent 
mechanical profiler with a higher resolution. Two hun- 
dred parallel profiles were recorded with 2050 points per 
profile (in which Ax = 32.5/j,m and Ay = 250/im). We 
assumed that the crack speed was constant which trans- 
lates in a linear relationship between position y and time 
t. As a consequence, we consider the one-dimensional 
profiles as descriptions of the advancing crack h(x,t). 
The complete spatio-temporal behavior of the surface can 
thus be obtained. 

In reference ||, the scaling form (||) was checked for 
the two first data sets. However, a careful inspection of 
the data collapse reported, Fig. 1 in ||, reveals that 
the scaling function goes like a power law for large ab- 
scissa values rather than be a constant. Also the slope for 
small values of the argument does not match well with 
the correct value. As we will see much better and more 
accurate results are obtained if, instead of assuming a 
Family- Vicsek behavior, we analyze the data on the ba- 
sis of an anomalous scaling. In this case, from (g) it is 
easy to see that the corresponding scaling function would 
be 

u-x ifu »i. < ( 9 ) 

in such a way that w(l,t)/l x — fl^OA )j where the 
label A denotes the anomalous scaling form. 
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In Figure 1 we present the data collapse of w(l,t)/l x 
vs. I ft 1 '* for the high resolution map of the crack sur- 
face obtained in the experiment. The best data collapse 
occurs for a global roughness exponent \ = 1.2 ± 0.1 
and z — 1.2 ±0.15. The same results were also found 
when the other two lower resolution data sets were used. 
The non-constant behavior for « C 1 in Fig. 1 is the 
main fingerprint of the anomalous character of the scal- 
ing. Figure 1 is in excellent agreement with a scaling 
function like (||). 

According to (^|) the power law u~ 0A1 for u <C 1 in Fig. 
1 corresponds to a local roughness exponent xioc = 0.79. 
This value can be compared with the value obtained from 
the long time behavior of the height difference correla- 
tion function, G(l,t) = ((h(x + /,£) — h(x, t)) 2 )^ 2 , that 
we plot in Figure 2 for the long times limit, t > P, and 
the highest resolution data. At long times, G(l) displays 
a power law behavior l Xloc that gives an independent de- 
termination of the local roughness exponent. As shown 
in Fig. 2, data fit to xioc — 0.79 in agreement with the 
value obtained from the slope —0.41 in Fig. 1. This esti- 
mate of xioc was confirmed by several other independent 
techniques: variable band width, return probability and 
wavelet analysis pT|]22| . 

We have also calculated the power spectrum and in 
Figure 3 we plot S(k, t) vs. k in a log-log plot for 
times t = 50,75 and 100 for the high resolution data. 
The curves are clearly shifted in time as corresponds 
to a power spectrum scaling function like Eq.(^) and 
not the Family- Vicsek one in (Q). S(k,t) decays with 
a power law fc -2 58 which is consistent with k~^ 2xioc+1 ^ 
and xioc = 0.79. We were unable to obtain a collapse of 
good quality form these data. 

Conclusions.- In this study we have shown that an 
anomalous dynamic scaling, (H) or ([?]), captures much 
better the features of the crack geometry than the stan- 
dard Family- Vicsek one. We obtained for an unsta- 
ble brittle fracture of a granite block a global expo- 
nent x = 1-2, which is different from the local exponent 
Xioc = 0.79. The very robust value of the local exponent 
has been demonstrated in the past. 

The existence of two different and independent rough- 
ness exponents that characterize the complete scaling of 
a surface, as it follows from Eq.(^), has important im- 
plications in the appearance (and geometry) of anoma- 
lously roughened surfaces. To illustrate this point we 
plot in Figure 4 two surfaces with exactly the same local 
j Xioc = 1/2, but a different global exponent, x — 3/4 and 
1/2. The two interfaces plotted in Fig. 4 are quite dif- 
ferent to the naked eye despite they have the same local 
roughness exponent. From Figure 4 it is clear that xioc 
necessarily gives just a part of the information about the 
scaling when anomalous roughening exists and indeed the 
global exponent x does give account of the large peaks 
taken by the surface in the anomalous case (i.e. when 
X Xioc)- The similarity with what happens to the 
patterns found in experimental cracks (see for instance 



Fig 1. in H) suggests that they may exhibit the same 
type of double scaling. More precisely, Fig. 1 in Rcf. U 
shows that the resulting surfaces in wood for tangential 
and radial fractures display strong different morphologies 
with the same local exponent \ioc 

~ 0.68. These authors 
found puzzling the large peaks taken by tangential frac- 
ture surfaces in clear contrast to the quite flat look of 
radial fractures, very much as occurs in the example we 
plot in Figure 4. Also Zhang et. al. p3[ in a numer- 
ical model of fracture in anisotropic materials found a 
local exponent roughly constant xioc — 0.7 surprisingly 
independent of the orientation of the fracture, although 
visible differences in the surfaces were noted (see Fig. 2 
in Ref. p3[]). Experiments in material showing differ- 
ent rupture modes by Bouchaud et. al. Q gave similar 
results. We believe that a determination of the global 
roughness exponent in all these experiments could allow 
a better characterization of the fracture surface morphol- 
ogy and its relationship with material properties. 

A second physical consequence of the anomalous scal- 
ing is that the saturation roughness is not only function 
of window size but also of the system size. Implications 
on the crack process are important. Information about 
the system size exists along the crack front during the 
propagation even at Rayleigh speed. This may illustrate 
the role of interactions between elastic waves and front 
propagation for the geometry of the front. 

More studies should be initiated to check weather or 
not the global roughness exponent is in general a valid 
index to characterize fracture surfaces. 
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and encouragement and R. Cuerno for a careful read- 
ing of the manuscript. JS thanks K.J. Mal0y for very 
fruitful discussions. This work has been supported by 
the MEC of the Spanish Government, the European 
Commission, the Ecodev-CNRS program and the GDR 
"Geomecanique des Roches Profondes" . 




loq^.t 



-2.0 -1.0 0.0 1.0 2.0 

log 10 ( l/t 1/z ) 

FIG. 1. Inset shows the local width vs. time for the higher 
resolution data of the crack experiment calculated over win- 
dows of sizes I = 100 (solid), 300 (dotted), 500 (dashed), 700 
(long-dashed) and 900 (dot-dashed) (all in units of the grid 
step Ax). Data are collapsed for \ — 1-2 and z = 1.2 in the 
main panel for windows sizes ranging from I = 10 to I = 1200. 
The non constant behavior at small values of the argument 
(main panel) reflects the anomalous character of the scaling, 
which agrees with Eq.(9). 
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FIG. 2. Height difference correlation function for the ex- 
perimental data displayed in Fig. 1. The straight line is a 
fit of the data and its slope 0.79 gives a determination of the 
local roughness exponent. 




FIG. 3. Power spectrum at times t = 50, 75, 100 calculated 
from the experimental data of higher resolution. The shift of 
the curves for different times is apparent and characteristic 
of an intrinsic anomalous scaling. The straight line has slope 
-2.58 and is in agreement with a the power decay k~^ 2xioc+1 ^ 
with xioc = 0.79. 



FIG. 4. Example of two fractal curves with the same local 
roughness exponent but different global one. One interface 
(solid) has \ = Xioc = 1/2 and is a realization of the simple 
dth = d^h + £ equation, where £ is a Gaussian white noise, 
and it is thus a true self-affine interface. The other curve 
(dashed) has \ = 3/4 and xioc = 1/2 and is a typical front 
of the random diffusion growth process [17,19], which is well 
known to exhibit anomalous roughening. 
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